ULM-TP/99-6 
HPL-BRIMS-1999-06 
September 1999 



Semiclassical form factor for chaotic systems with 
spin 1/2 



Jens Boltefflj] and Stefan Keppeler^ 

f Abteilung Theoretische Physik, Universitat Ulm, Albert-Einstein- Allee 11, D-89069 
Ulm, Germany 

\ School of Mathematics, University of Bristol, University Walk, Bristol BS8 1TW, 
United Kingdom 

§ BRIMS, Hewlett-Packard Laboratories, Filton Road, Stoke Gifford, Bristol 
BS34 8QZ, United Kingdom 

Abstract. We study the properties of the two-point spectral form factor for 
classically chaotic systems with spin 1/2 in the semiclassical limit, with a suitable 
semiclassical trace formula as our principal tool. To this end we introduce a regularized 
form factor and discuss the limit in which the so-called diagonal approximation can be 
recovered. The incorporation of the spin contribution to the trace formula requires an 
appropriate variant of the cquidistribution principle of long periodic orbits as well as 
the notion of a skew product of the classical translational and spin dynamics. Provided 
this skew product is mixing, we show that generically the diagonal approximation of 
the form factor coincides with the respective predictions from random matrix theory. 



PACS numbers: 03.65. Sq, 05.45.Mt 



|| E-mail address: bol@physik.uni-ulm.de 
% E-mail address: kep@physik.uni-ulm.de 

+ Address after 1 October 1999: Abteilung Theoretische Physik, Universitat Ulm, Albert-Einstein- 
Allee 11, D-89069 Ulm, Germany 



Semiclassical form factor for chaotic systems with spin 1/2 



2 



1. Introduction 

One of the major paradigms of quantum chaos is the conjecture of Bohigas, Giannoni 
and Schmit (BGS) [Tj which states that the local statistics of energy spectra of 
(generic) individual quantum systems, whose classical analogues exhibit (strongly) 
chaotic behaviour, can be well described by that of ensembles of large random matrices. 
The symmetry properties of the relevant matrix ensembles have to be chosen according 
to the symmetries of the quantum system under consideration. In case the system 
is invariant under time reversal and has integer total angular momentum its local 
eigenvalue statistics are conjectured to be that of the Gaussian Orthogonal Ensemble 
(GOE). If time reversal invariance is broken one expects local statistics according to the 
Gaussian Unitary Ensemble (GUE). 

However, if the total angular momentum of the system is half-integer and the system 
is invariant under time reversal all eigenvalues show Kramers' degeneracy 0, [| and their 
statistics have to be compared with the Gaussian Symplectic Ensemble (GSE). In the 
GOE- and in the GUE-case there is plenty of numerical evidence available in favour of 
the BGS-conjecture, see, e.g., H3» S 0' whereas only few examples have been studied in 
the GSE-case as, e.g., in || (7[. For quantum systems whose classical limit is integrable, 
i.e., which shows regular behaviour, one expects the local eigenvalue statistics to follow 
the laws of a Poisson process M. For the analytical treatment semiclassical methods, 
in particular semiclassical trace formulae, have become the most important tools since 
Berry and Tabor || investigated the behaviour of the spectral form factor for classically 
integrable systems by means of an appropriate trace formula. By making use of the 
Gutzwiller trace formula ||, [|, |lTj, [TT], [12| , Berry provided a semiclassical theory for the 



spectral form factor fl3| of classically chaotic systems without spin. Based on the so- 
called diagonal approximation he could explain the semiclassical asymptotics of the form 
factor for small values of its argument, thus recovering the GOE- and GUE-behaviour, 
respectively. 

In this paper our aim is to show that Berry's semiclassical treatment of the two- 
point form factor can be carried over to quantum systems with spin 1/2, whose classical 
translational dynamics are chaotic. We base our analysis on the semiclassical trace 
formula for the Dirac equation that we developed recently [TJ], In this trace formula 
the presence of spin is reflected in a modification of the amplitudes with which the 
periodic orbits of the translational dynamics contribute. This modification arises from 
a spin dynamics that involves a 'classical' spin precessing along the periodic orbits. 
The central part of this paper therefore consists of calculating the effect of this spin 
contribution to the semiclassical form factor. For our analysis we use a two-point form 
factor whose definition differs slightly from the one that is more commonly used in 



spectral statistics as, e.g., in fI3|. We rather prefer the point of view adopted in [TBI, [TT 



Both definitions, however, are equivalent in the limit where infinitely many eigenvalues 
are taken into account. We also stress that both the form factor and the associated 
correlation function are distributions and hence have to be evaluated on suitable test 
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functions. This approach makes a spectral average obsolete and enables one to state 
the BGS-conjecture, specialized to the form factor, in a precise manner. Moreover, the 



lacking self- averaging property discussed in [18] poses no difficulty in this context. 

The paper is organized as follows. In section ^ we introduce our definition of the 
spectral form factor for a finite part of the spectrum and discuss in which sense one 
can expect to recover the form factors given by random matrix theory. Section |3] is 
devoted to the definition of a regularized form factor which we evaluate semiclassically 
using trace formulae for the Dirac as well as for the Pauli equation |14|, [TI|. We also 
invoke the diagonal approximation and briefly discuss its range of validity. In section f| a 
suitable version of the equidistribution principle of long periodic orbits is used in order 
to obtain the semiclassical asymptotics of the diagonal form factor. In this context 
we employ the notion of a skew product of the translational and the spin dynamics, 
the ergodic properties of which determine the semiclassical asymptotics. Our principal 
results are then summarized in section 5. Namely, depending on the presence or absence 
of quantum mechanical time reversal invariance, and provided the dynamics of the skew 
product is mixing, we can recover a GSE- or GUE-behaviour of the diagonal form factor, 
respectively. The relation between classical and quantum mechanical time reversal as 
well as the equidistribution of long periodic orbits are discussed in two appendices. 

2. The form factor for quantum systems with half-integer spin 

The two-point correlations of a discrete quantum spectrum are conveniently measured 
by either the two-point correlation function R 2 or by the two-point form factor K 2 , 
which is related to R 2 through a Fourier transform. Before defining these quantities 
one usually unfolds the spectrum, i.e., the eigenvalues Ej. are rescaled to x k such that 
the unfolded eigenvalues have a mean separation of one. This means that the spectral 
density d(x) of the unfolded spectrum allows for a separation 

d(x) := ^6(x-x k ) = l + d fl (x) , (2.1) 



px+Ax j 

/ d fl (y)dy=——#{k;x-Ax<x k <x + Ax}-l (2.2) 

Jx-Ax ILAX 



such that 

px+Ax 

2Ax 

vanishes as x — > 00, Ax — > 00, Ax/x —> 0. For a finite part of the spectrum, containing 
N unfolded eigenvalues enumerated as xi, . . . ,xn, one defines the two-level correlation 
function by 

k,l<N 

Accordingly, the two-level form factor is defined as 

K 2 (t; N) := I R 2 (s; N) e" 2 " iTS ds = e~ 2nW{xk - Xl) - 8{t) . (2.4) 
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Since both quantities are distributions, which is most clearly seen in the case of the 
correlation function (12. 31), one should evaluate these on test functions d> G 



/ K 2 (t;N)<P(t) dr = -J- V f 0(r) e^^^dr - 0(0) 

N k,l<N J * 

-j- ^ 0(2tt(^ - a*)) - 0(0) (2.5) 



k,KN 



R 2 {s;N)(f){2irs) ds . 

We remark that since the form factor is obviously even in r, it suffices to consider only 
even test functions 0. The convention for the Fourier transform that was used, and that 
will be used in all of what follows, is 

/(*) = / f{x) e ixk dx and f(x) = -L / f(k) e"^ dA; . (2.6) 
Jr ^ Jr 



After smearing with a test function all expressions occurring in Q2.5Q are obviously finite. 
In this form a semiclassical analysis of either the form factor or the correlation function 
can be carried out. If is chosen non-negative, the left-hand side of fl2.5|) can also 
be viewed as the mean value of K 2 (t; N) when r is drawn randomly with probability 
density 0. Since we always understand the form factor in the above sense, the absence of 



a self- averaging property discussed in |L8[ is not essential for our further considerations. 

For a given quantum Hamiltonian H the unfolding of its discrete spectrum shall 
proceed in the following manner. We consider a spectral interval 

/ = I(E, h) := [E — fajj, E + fajj] , u > , (2.7) 

that has no overlap with a possible essential spectrum of H. Then the condition Ek G / 
is equivalent to 

Ei. — E , . 

-lu<— < uo . (2.8) 

n 

The number of eigenvalues contained in J, 

Nj := E k G 1} , (2.9) 
can be estimated semiclassically as 

iVj ~ 2hLud(E) = -T H (E) , ft -> , (2.10) 

7T 

where denotes an appropriate mean spectral density and Tjj(E) := 2irhd(E) is the 
Heisenberg-time. A convenient definition of d(E) can be derived from the semiclassical 
trace formula for H in that it shall denote the contribution coming from the singularity 
of Tr exp[— j^Ht] at t — to all polynomial orders in h see, e.g., |T^, |T1], |T2|, |T5], [T9|. The 



spectra that we are going to consider below are such that in the semiclassical limit h~ — ► 
the Heisenberg-time tends to infinity, Th — > oo. For example, given in the gap of the 
essential spectrum of a Dirac-Hamiltonian Hp, i.e., in typical cases —mc 2 < E < mc 2 , 



the mean spectral density reads [14, 15 



7, „s _ VOl Qt + VOl Q TP r _ „,,v, ,_..s 
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where Q E denote the hypersurfaces of energy E in phase space corresponding to the 
classical Hamiltonians 



if^p, x) = &p(x) ± y (cp - eA(x)f + m 2 c 4 (2.12) 

for relativistic particles of positive and negative kinetic energy, respectively, in the static 
external electromagnetic fields generated by the potentials ip and A. Therefore, in the 
semiclassical limit the number of eigenvalues in the interval / increases, although its 
length |/| = 2fkv shrinks to zero. A completely analogous argument applies to Pauli- 



Hamiltonians 15 



We now define the unfolded spectrum through 

x k := E k d{E) and x := Ed(E) . (2.13) 

The condition E k G I(E,h) is hence equivalent to x k G [x — Ax,x + Ax], where 
Ax := ^Th- With this choice indeed x — > oo, Ax — > oo, such that Ax/x = frw/E — > 
in the semiclassical limit. In this context the quantity ( |2.2|) reads 



— - #{k; x-Ax<x k <x + Ax}-l = ^-Nj-l, (2.14) 
2Ax oj±h 



such that (|2.10|) ensures its vanishing in the semiclassical limit. 

From now on we will choose the numbering of the eigenvalues E k and x k , 
respectively, in such a way that the eigenvalues in I are given by 

E 1 <E 2 <...<E Nl . (2.15) 

Changing the value of h therefore alters the numbering of the eigenvalues. As a 
consequence, the condition E k G I is equivalent to k < Nj. At this place we recall 
that the semiclassical limit H — > 0, or Tjj — > oo, implies Nj — > oo. For the form factor 
we now obtain 



2 (2.16) 



iVj 



Ei E k — E 
X[ 



Nj ^ ^ l ' J V h 

k 



h 



'-rT H E k 



-St) 



where the A;-sum extends over all eigenvalues of H and X[-u,u] denotes the characteristic 
function of the interval [—a;, a;] that occurs due to the condition (|2.8|). 

In the preceding discussion we tacitly assumed that the discrete spectrum of 
the quantum Hamiltonian carries no systematic degeneracies. In order to achieve 
such a situation one has to remove all symmetries. As opposed to geometric or 
internal symmetries that have to be realized by unitary representations of the respective 
symmetry groups, the time reversal operation must be implemented by an anti-unitary 
operator T, see [[§ and Appendix A| . For single particles of spin s the square of T 
depends on s being integer or half-integer in that T 2 = (— l) 2s . In case the quantum 
system is time reversal invariant, i.e., [H, T] = 0, and has half-integer spin this leads to 
Kramers' degeneracy [0: Since T 2 = —1 implies that every vector ip ^ in the Hilbert 
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space is orthogonal to Tip, all eigenvalues of H are (at least) two-fold degenerate, see 



2lf for details. 

Following the usual practice, we will remove Kramers' degeneracy by replacing each 
degenerate pair E 2k = E 2 k+i of eigenvalues by one of its representatives. Thus the mean 



spectral density is lowered by a factor of two. In analogy to Q2.13| ) the unfolding i— > x\. 



of the so modified spectrum can therefore be achieved through the choice x k := 
The modified form factor then reads 

K 2 (r;N) = — e" 2 " iTfe - £i) - 5(t) 

k,l<N 
k,l odd 

= J_ V e- w i( Jrt ) --S(-\ ( 2 - 17 ) 
2N ^ 2 \2) 

k,l<N 

In this setting the conjecture of Bohigas, Giannoni and Schmit |1[ states that for 
individual (generic) classically chaotic quantum systems of particles with half-integer 
total spin and no unitary symmetries one should obtain 

lim [ K 2 {T;N)(f){T) dr = f K^ SE {t) ^{t) dr (2.18) 

for all test functions <fi £ Cq°(]R). Here K 2 SE denotes the two-point form factor of the 
Gaussian symplectic ensemble of random matrix theory, 

^otp I hW\ — t|t| losr 11 — Irll for Irl < 2 , 

K% SE (r) = < 1 " II-' / 2 _ 1Q 

2 1 for Irl > 2 , 



see, e.g., pi|. When time reversal invariance is lacking, the respective conjecture reads 



! 



lim / K 2 (t; N) <j)(r) dr = / fQ^(r)0(r) dr , (2.20) 

TV— >-oo 



where now the form factor of the Gaussian unitary ensemble [21] should occur 



k? ve {t) = r ' 1 7 (2.21) 





for r 


[1 TI 


for r 



In our subsequent semiclassical investigations we will in both cases, i.e., with and without 
time reversal invariance, consider the form factor K 2 (t; Ni) as it is given in ( |2.16| ). When 
dealing with the case of time reversal invariance we appeal to the relation fl2.17|) . 



3. The semiclassical form factor 

Since the work of Berry and Tabor || on the distribution of eigenvalues for classically 
integrable systems, semiclassical trace formulae have found numerous and fruitful 
applications in the analysis of spectral statistics. A prominent example is Berry's 



analysis of the spectral rigidity [13|, which relies in an essential way on a semiclassical 



evaluation of the two-point form factor based on the Gutzwiller trace formula. In this 
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work it already became apparent that present semiclassical methods at most allow to 
study the form factor in the restricted range \r\ < 1, see also []I9[ for a review. Only 
recently, improved techniques have been developed that might allow to extend the 
semiclassical analysis of spectral statistics. In this paper, however, we follow the more 
traditional path in that in the end we consider the so-called diagonal approximation for 
the form factor. 

The two-point form factor as given in ( [2.16D requires to establish a trace formula 
for the sum 

(3-D 

However, the general structure of (convergent) semiclassical trace formulae, see, e.g., 
|10| , 11, 12, 15, 19], necessitates the use of a smooth test function p £ C°°(R) with 
Fourier transform p £ C^°(M). One therefore has to replace the sharp cut-off, provided 
by the characteristic function in ( |3.1| ), by a smoothened substitute. For this reason we 
now introduce the regularized form factor 

2 

-S(t), (3.2) 



k 



where rj £ C°°(R) is a test function with Fourier transform fj £ C£°(R), but that 
is otherwise arbitrary at the moment. Later we will introduce a further normalization 
condition. In the following we will consider both relativistic and non-relativistic particles 
with spin 1/2. In the relativistic case, when one is dealing with a Dirac-Hamiltonian Hp, 
the function x £ Co°(R)> which is not to be confused with the characteristic function 
X[-u>,uj\i is necessary to truncate the essential spectrum of Ho- In typical situations x 
should therefore be supported in the interval (—mc 2 ,mc 2 ), where the eigenvalues of 
Hd are located. When these do not accumulate at some point, one could also leave out 
the truncation x from ( |3.2| ). 

We are now in a position to use the test function 

p{e):=7 1 {e)e- 2 ^ E ^ +E ^ (3.3) 

in the semiclassical trace formula for the Dirac equation that was developed in |l4l [if 



h \ 



rp (3.4) 

7 k^O 

The outer sum on the right-hand side extends over all primitive periodic orbits 7 of 
energy E, with periods T 7 , of the two classical flows generated by the Hamiltonians 
( |2.12j ). The inner sum then is over all fc-fold repetitions of primitive orbits, formally 
including negative ones. The weight attached to each pair (7, k) reads 

A, k := ^4 r ei^W-qk^ m + Q/^i / 3 5 x 

7 ' |det(M*-l)|3 1 ^ 
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Here c? 7 G SU(2) denotes the semiclassical time evolution operator for the spin degrees of 
freedom along the primitive periodic orbit 7 of the translational dynamics. Furthermore, 
S^(E) denotes the action of 7 and /i 7 is its Maslov index. The (monodromy) matrix M 7 is 
the linearized Poincare map transversal to 7. In the form given in ( |3.4j ) the trace formula 
is valid for all cases where the classical flows have only isolated and non-degenerate 
periodic orbits. An analogous trace formula, with appropriate simplifications, is also 
available for Pauli-Hamiltonians, see fll5| . 

Upon choosing the test function (fDf ) in the trace formula (|3.4|) , its left-hand side 
reads 

£ x{E t ) P (^S) = £ x(E k ) V (^) e~i™> , (3.6) 
k ^ ' k ^ ' 

and is hence the appropriate starting point for a semiclassical analysis of the form 

factor, compare (|3.2j ). As a first ingredient on the right-hand side of ( |3.4|) one requires 

the Fourier transform of the test function (|3.3|) , which is given by 

p(t)=e- L n ETT "fi(t-rT H ) . (3.7) 

For convenience we now choose rj to be even and real-valued, which implies that 17 
also shares these properties. Furthermore, the truncation \ of the essential spectrum 
shall be such that x[E) = 1. Thus, the trace formula yields the following semiclassical 
representation of the regularized form factor, 

K™(r; T H ) = - 5(r) + [f,(rT H )} 2 [1 + 0(h)] 

+ E E ^ ^(tTh) v(kT, - tT h ) A^ k [1 + O(^)] (3.8) 

+ ^- S~] y~] Till. fj^kTry - tT h ) t)(tT h - k'Ty) A lik A 7 > k , . 

" 7,7' k,k'y£0 

In a next step we are going to test the semiclassical form factor with some <fi G Cg°(R), 
compare (21S). To this end one needs the integral 



F(T, T) := [ <f>(r) f)(T - tT h ) fj(rT H - T) dr 
Jr 



(3.9) 



whose leading term in the semiclassical limit T# — > 00 can be calculated by introducing 
the Fourier representations for and fj. A straight-forward calculation then yields 



F(T, T) = ±- [ [ [ 4>(t) r)(e) v(s') e^ T ~^ sf^ + s-s') de' de dt . (3.10) 
J-h JrJrJr \ + h J 



Changing variables from e,e f to u := e' — e and v := (e f + e)/2, and employing the 
expansions 

■„U^ ; )=, l (v) + o(±) , (3.11) 
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finally shows that 

F(T, T>) = ±- I [ }(t) e- it( ^ ^vf dvdt + O ( * 

J-h JrJr \ j h 

r)*r)(T-r') + 0^ , 



1 

7^ 



T + T r 



(3.12) 



where the convolution 

fj * r)(t) := / f)(t - t') fj(t') dt' = 2n I 7](vf e ltv dv 
Jr Jr 

enters. We therefore conclude that 



(3.13) 



Kf\r-T H )ct>(T) dr = m 



-1 + 



Atxuo 



fj * 77(0) 



+ o(h) 



7 fc^o 



(3.14) 



7,7' k,k'jtO 
X 



47ru; 



7' 



kT-^i ~\~ fcj Try! 

27^ 



A-7,fc -A7',— fc' 



At this point we introduce the normalization of 77 announced previously. Guided by the 
simple observation 



(3.15) 



we require the same normalization for the smooth substitute 77 of the sharp cut-off 

A.[— U>,tJ\ ) 



— fj * ?)(0) = / 77(e) 2 de = 2u . 



(3.16) 



As a consequence, the leading semiclassical order of the first line on the right-hand side 
of ( |3.14| ) vanishes. Furthermore, since the Fourier transform fj of the test function rj 
is required to be compactly supported, the second line is a finite sum, multiplied by 
1/Th- A similar argument applies to the third and fourth line, apart from the diagonal 
contribution with /cT 7 = k'Ty , where f) * ry(0) occurs and thus no such cut-off is present. 
Due to the above reasoning it is tempting to assume that in the semiclassical limit 



Th — > 00 the right-hand side of ( |3.14j ) is completely fixed by the contribution of the 
diagonal form factor 



J- T-r 



7,*| 2 <M^-^) • (3-17) 

H 7 fc^O 

Here we assumed that k'Ty = kT^ implies Ay^i = A^, see [Appendix A for the spin 
contribution, and then g^^ denotes the number of pairs (7, k) such that /cT 7 has a given 
value. It is, however, well known that the above assumption is not justified. The reason 
for this lies in the subtleties of the limits involved. In order to arrive at the left-hand 
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side of ( [2.18 ), or of ( 2.20|) , one first has to remove the smoothening of the characteristic 



function X[-u>M in that a sequence of functions n G C°°(E) approaching X[-w,w] has to be 
considered. Then, in the limit, 17*17 is no longer compactly supported. Indeed, according 
to ( |3.13| ) one obtains fj * fj(t) = 4n sm(ut)/t. Still, the periodic-orbit sums in ( |3.14| ) 



are truncated by the test function 0. However, in the semiclassical limit this cut-off 
is being removed. Moreover, for long periodic orbits the differences kT 1 — k'Ty can 
become arbitrarily small so that fj * fi{kT 1 — k'Ty) only provides a modest truncation of 
near-diagonal contributions. The only example where it could be rigorously shown [[U| 



that the diagonal form factor itself produces the correct limit, if the test functions are 
restricted to those that are supported in the interval [—1, 1], is that of the correlations 
of the non-trivial zeros of principal L-functions, including the case of the Riemann zeta 



function, see also p3| . Rudnick and Sarnak [|16|] even proved an analogous result for 
general n-point correlations. 

As announced previously, we now invoke the diagonal approximation, i.e., we leave 
aside the contribution of /fg'^Cr; T#) — K^ 9 {t\Th) to ( |3.14|) . This procedure, which 
goes back to Berry [|13j, is generally supposed to reveal the correct behaviour of the 



form factor for small |r|. The reason for this being that if the test function is 
supported in a small interval, the contribution of long periodic orbits to (|3.14| ) is 



truncated. Furthermore, due to the normalization ( |3.16| ) the diagonal form factor is 



independent of the smoothening rj. This convenient fact exempts one from the need to 
discuss the removal of this smoothening. In order to test now the range of small |r| 
one should restrict the class of test functions to those supported in intervals [— r', r'], 
where r' > is 'small enough'. Recalling that K^ a9 and are even in r an integration 
by parts yields 



Kf a \r- T H ) 0(r) dr = 2 / #f aff (r; T H ) 0(r) dr 



o 



2 (3-18) 

1 H 



What is required now is the asymptotic behaviour of the periodic-orbit sum in (|3.18| ) 
as Tjj — > oo. Since the contributions of repetitions of primitive periodic orbits are 
asymptotically suppressed due to their stronger instabilities, compare also ( fOD below, 
in the following we only take the k = 1-term of the sum over the repetitions into account. 
As a consequence we therefore have to study the asymptotics of the periodic-orbit sum 

|det(M 7 -l)| 1 ] 

in the double limit Th — > oo, r — > 0, such that tTh — > oo. In order to simplify this 
task we now make two assumptions, which should be verified in all cases that could be 
considered as 'generic' in any reasonable sense: 

(i) The periods T 7 of primitive periodic orbits shall be such that any finite subset of 
them is linearly independent over Q. This implies that the multiplicities g 7i fc are 
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independent of k, i.e., g 1 ^ = g 7 - 

The subset of primitive periodic orbits 7 with g 1 7^ g is of density zero in the set 
of all primitive periodic orbits, 

where g = 2 in case the classical dynamics are time reversal invariant, and g = 1 
when time reversal symmetry is absent. In this context time reversal invariance 
does not only mean that an orbit 7 is geometrically identical to its time reversed 
partner, but also that both orbits yield the same contribution of the spin degrees of 
freedom to the trace formula, which then implies that is invariant under time 
reversal. In [Appendix A| we show that this condition is a consequence of quantum 



mechanical time reversal invariance. 

Under these assumptions the factors g 7i £ can be replaced by g and can then be pulled 
out of the sum fl3.19|) . 

4. Classical periodic-orbit sums and the contribution of spin 

The aim of this section is to obtain the leading semiclassical behaviour of the periodic- 
orbit sum (|3.19|) . Apart from the appearance of the Heisenberg-time only quantities 
related to the classical flow enter this expression. It therefore seems appropriate to 
invoke results about the distribution of periodic orbits in phase space. In order to retain 
a certain convenient generality, we will not specify the classical translational dynamics 
further, except for the following assumptions: 

(i) The classical flow & H : fig -» flg on the compact hypersurface Qe of energy E in 
the 2(i-dimensional phase space is generated by some Hamiltonian function H(p, x) 
and hence preserves the (normalized) Liouville measure 

dfi E (p,x) := — j-pp 5(H(p, x) — E) d d pd d x (4.1) 
vol He 

on Q E - 

(ii) & H is ergodic with respect to Liouville measure. 

(iii) is hyperbolic on all of VLe- 

In the case of the semiclassical form factor for a Dirac-Hamiltonian these requirements 
shall apply to both classical flows, i.e., to those generated by the two classical 
Hamiltonians if ± given in ( gT2|) . Moreover, we now assume that for a given energy E 



there will only be either a contribution coming from the dynamics generated by H + or 
from the dynamics generated by H~ , but never from both at the same time. This is not 
a strong restriction since it only excludes situations in which Klein's paradox |24 can 
appear. 

The hyperbolicity of the classical flows implies that in particular all periodic orbits 
are either hyperbolic or loxodromic. This means that all monodromy matrices M 7 have 
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eigenvalues with moduli strictly different from one. Since the eigenvalues occur in pairs 
of mutually inverse numbers, we denote them as e ± ( u ~<- j+w ~>' j \ u 7i j > 0, v y j G [0, 2ir), 
j = 1, . . . ,d — 1. Thus 



det(M 7 - 1)1 = IJ I (e" 7 '^^ - l) (e " ir - l) | 




nil 

j=l / j=l 



d-l 

2 



(4.2) 



U J,3 



The stability exponents u Jt j are related to the Lyapunov exponents A 7 j of 7 through 
u 7 j = A 7J T 7 so that u y j — > 00 as T 7 — > 00. Hence, in this limit one obtains 

p 7 := exp f J • (4.3) 



det(M 7 - 1)| 



Since the semiclassical limit of the periodic-orbit sum ( |3.19| ) is dominated by the 
contribution of long periodic orbits, ( [4.3[ ) allows to analyze ( ^.19| ) in terms of periodic- 
orbits sums that are familiar from equidistribution theorems of periodic orbits, see, e.g., 



25 



For the kind of Hamiltonian flows characterized above mean values of observables 
on Qe with respect to Liouville measure can be calculated with the help of appropriate 
periodic-orbit sums. We postpone a detailed discussion of this matter to |Appendix B| , 



from which we here only quote that for any continuous observable a one obtains the 
representation 



:= j a(p,x) dfi E (p, x) = lim - ^ T 7 a 7 p 7 , (4.4) 



n E 7,T 7 <T 

where a 7 denotes an average of a along the periodic orbit 7, 



a 7 :-- 



1 f 7 

— / a ($* H (p,cc)) dt , with (p,x) E 7 • (4.5) 

J 7 JO 

We remark that a heuristic derivation of an analogous identity to (4.4) was given by 
Hannay and Ozorio de Almeida p6| . 

One can apply the relation ( |4.4|) to determine the leading semiclassical behaviour of 
( P-19| ) once one has chosen a suitable observable a whose average a 7 yields the quantity 



(tr<i 7 ) 2 appearing in fl3.19| ). This, however, can only be achieved in an indirect manner. 



Our choice of the observable requires to recall the semiclassical time evolution of the 
spin degrees of freedom along the trajectories of the classical flow Let d(p, x, t) G 
SU(2) denote the solution of the spin transport equation |T4] , |T5] 

d(p,x,t)+iM(& H (p,x)) d(p,x,t) = , d{p, x, 0) = t 2 , (4.6) 

where the time derivative is understood along the trajectory $^(p, x). M is a certain 
hermitian and traceless 2 x 2-matrix valued function on Qe, whose precise form depends 



on the quantum Hamiltonian under consideration, see |L5| for details. Geometrically. 



d{p, x, t) can also be interpreted as a parallel transporter in some vector bundle so that 
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d(p, x, T 7 ), with (p, x) G 7, is the holonomy associated with the periodic orbit 7. Since 
its trace is invariant under a shift of the initial point (p, x) G 7 along the orbit, one can 
introduce the notation tr<i 7 := tr d(p, x, T 7 ). We are thus in a position to define the 
observable 

a(p,x,t) := [tr d(p,x,t)} 2 , (4.7) 

which is a function on fi^ that in addition depends on a parameter t. Due to the above 
remark concerning the interpretation of d(p,x,T^) as a holonomy, the average ( [4.5|) of 
this observable along a periodic orbit 7, when t = T 7 is chosen, yields 

a^(T 7 ) = yJ {faiP^T^) dt' = (trd 7 ) 2 , (4.8) 

for any (p, x) G 7. Without the choice t = T 7 , however, a 7 (t) is not related to (tr<i 7 ) 2 . 
We can hence now employ ( |4.4| ) to deduce the asymptotic relation 

a\t) Vl ~ It 2 a s (t) , T - 00 , (4.9) 

7,T 7 <T 

which is valid for any t. Notice that here we have introduced an extra power of T 7 in 
the same manner as in (|B.7|) -( |BT8l) . We now differentiate with respect to T, 

^T 7 2 a 7 (t)p 7 5(T-T 7 )~Ta E (t), T -> 00 , (4.10) 

7 

and then choose t = T. Together with ( [4.8| ) this allows to conclude that 

^T 2 (tr rf 7 ) 2 Vl 8(T - T 7 ) ~Ta E (T) , T -> 00 . (4.11) 

7 

Thus, at this point we have obtained the asymptotic relation 

K^^T^^gra^rTH) (4.12) 

for the diagonal form factor ( |3.17|) in the limit Tjj —>■ 00, r — ► such that tTh — > 00. 

The remaining task therefore consists of determining the asymptotics of a E (T) as 
T — > 00. In order to achieve this one has to go back to the representation 

a E {T)= [ [ti d(p, x,T)} 2 dfx E (p,x) (4.13) 

of a E (T) as an average over phase space. Since the T-dependence involves d(p, x, T) one 
might anticipate that the limit T — > 00 of ( |4.13|) depends on certain ergodic properties 



of the spin dynamics. The latter being considered along trajectories of the translational 
dynamics, one hence has to combine both dynamics in a suitable way. In ergodic 
theory the relevant construction is known as a skew product, see, e.g., ||27|| . Appropriate 
ergodic properties of the skew product dynamics will then allow for a determination of 
the asymptotic behaviour of (|4.13|) . Let us therefore now construct the skew product 
of translational and spin dynamics. To this end one defines a flow Y* on the product 
phase space M. := VLe x SU(2) in the following way, 

F'((p, *),<?) := (& H (p,x),d(p,x,t)g) (4.14) 
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for (p,x) G Vl E and g G SU(2). The initial condition Y°((p,x), g) = ((p,x),g) is 
obviously fulfilled, and the composition law Y t+t ' = Y l o Y 1 ' immediately follows from 
the relation 

d{p, x, t + t') = d(& H (p, x),t) d(p, x, t') (4.15) 

that can be concluded from (|4.6|) . On M. one then defines the direct product fi := 
He x /J-h of Liouville measure \ie and of the normalized Haar measure //# of SU(2). We 
recall that the latter is the unique normalized left- and right-invariant positive Radon 
measure on the group manifold, see, e.g., HI. Due to both the invariance of Liouville 
measure under the Hamiltonian flow $^ and the left-invariance of Haar measure, the 
product measure /i is invariant under Y l . A dynamical system of this kind is known as an 
SU(2)-extension of $^ or, more generally, a skew product. The spin dynamics defined 
by (|4.6|) is then called a cocycle for $^ with values in SU(2). For further information 
see, e.g., |27]. 

In addition to the assumptions made for & H in section |3], in the following we will 
assume that Y l is (strongly) mixing. This implies that for any F G L 2 (A4 x Ai) 

lim / F(Y t ({p,x),g),((p,x),g)) dfi((p,x),g) 

^ (4.16) 
= / / F{((p,x),g),((Z,y),h)) d t i({p,x),g)d^,y),h) . 

JM JM 

We remark that usually the mixing property is defined for pairs of functions F\,F% G 
L 2 (j\4). However, if one views L 2 (M x M) as L 2 (M) ®L 2 (M) and introduces a tensor- 
product basis, (|4.16 ) follows immediately because every element of this basis fulfills the 
usual mixing property. If now the function F does not depend on the translational 
degrees of freedom, i.e., F : SU(2) x SU(2) — > M, the mixing property (|4. 16|) yields 

lim/ F{d{p,x,t)g,g) d»{(p,x),g) = / / F(g, h) dfx H (g) dfx H (h) . (4.17) 

^^JM JSU(2) JSU(2) 

A suitable choice of the function F then allows to determine the asymptotic behaviour 
of d E (T) as T —>■ oo from ( |4.17| ). In order to achieve this we recall the representation 
( |4.13| ) of d E (T), which obviously can also be written as 

a E (T) = [ [ [tr(d(p, x, T)gg- 1 )] 2 dfj, E (p, x) d^ H {g) . (4.18) 
isu(2) Jn E 

Upon now defining the function F(g,h) := [tr^/i -1 )] 2 , g,h G SU(2), one can employ 
Ql.l7| ) to conclude that 

lim a E (T)= lim / F (d(p, x, T)g, g) dfx{{p,x),g) 

JM (4.19) 



Mgh- 1 )} 2 d^ H {g) d/j H (h) . 

'SU(2) JSU(2) 

Substituting g' = gh~ x in the inner integral and using the right-invariance of the 
integrand does not depend on h any more. Thus we obtain 

lim a E (T) = [ [trgf dfi H (g') , (4.20) 

T ^°° JSU(2) 
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i.e., in the limit T — > oo the expectation value of [tr d(p, x, T)] 2 , when averaged over 
phase space, can be computed by an average over the group SU(2) with respect to Haar 
measure. 

The same holds obviously true for any moment of tr d(p,x,T) so that the 
asymptotic distribution of tr d(p, x, T), when the initial points (p, s) 6 are uniformly 
distributed with respect to fiE, can be computed via 

lim he {{p, x) G Qe] trd(p, x, T) G [a, b}} 

T — >oo 

rb r (4.21) 

/ 6(tvg-w) d[M H (g) dw . 

'a JSV(2) 

In order to evaluate the integral over SU(2) explicitly we remark that any g G SU(2) 
can be represented as 

3 

g(u) = u Q ±2 + icrn with u = (uo, u) G 1R 4 and V^tt| = l. (4.22) 

3=0 

In this parameterization the Haar measure is given by (see, e.g., [p8 |) 

d» H (g(u)) = 1 8 ^ - 1 j d ' u ■ ( 4 - 23 ) 

A simple calculation now shows that the distribution (|4.21|) obeys a semicircle law, i.e., 
its density reads 



p(w)= I S(trg-w) dfA H (g)= { n V (2) 2 < W < +2 , ^ 



SU(2) 



else 



What is required in (|4.20|) is the second moment of the distribution (|4.21 ). With the 
help of (|4.24|) this can now easily be computed to yield one, i.e., 

d E (T) ~ 1 as T -> 00 . (4.25) 



The integrated version of ( 4.11 ) therefore reads 

J2 T 7 2 (trcg 2 p 7 ~V , T^oo. (4.26) 

7,T 7 <T 

Furthermore, ( |4. 12j) and ( |4.25|) imply the asymptotic behaviour 

Kt 9 (r;T H )~g T (4.27) 
of the diagonal form factor in the regime T# — > 00, r — > 0, rT^ — > 00. 



5. Summary and Conclusions 



In section [3] we demonstrated how the two-point form factor for a quantum system with 
spin 1/2 can be analyzed semiclassically by making use of an appropriate trace formula. 
The diagonal approximation led us to the periodic-orbit sum (|3.19|) whose asymptotics 
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determines the behaviour of the diagonal form factor for small r in the semiclassical 
limit. Our conclusion (|4.27|) now has to be compared with the relevant random matrix 
results ( |2.19 ) and ( |2.21| ). To this end we have to distinguish time reversal invariant 
quantum Hamiltonians from non-invariant ones. But let us first summarize our findings. 
In |Appendix A| we show that quantum mechanical time reversal invariance implies that 
not only the classical translational dynamics are time reversal invariant but, moreover, 
the spin dynamics behave in such a way that also the amplitudes appearing in 
the semiclassical trace formula are invariant under time reversal. This leads to the 
occurrence of the multiplicities g^^ m the expression (|3.19|) . Then, when the generic 
conditions stated at the end of section ||] are met, we can pull out the factors g^^ from 
the sum and replace them by either g — 2, if (quantum mechanical) time reversal 
invariance is present, or else by g = 1. If, furthermore, the translational dynamics are 
ergodic and hyperbolic, the equidistribution of periodic orbits allowed us to determine 
the spin contribution to the amplitudes A ljk . We further requested the skew product of 
the translational and the spin dynamics to be mixing. This then enabled us to identify 
the distribution of the traces of the spin-transport matrices and to calculate its second 
moment, which enters through the amplitudes Ay^. 

Hence, if quantum mechanical time reversal invariance is absent, the diagonal 
approximation ( [4.27] ) states that K^ 9 (t]Th) ~ t, which is identical with the diagonal 
approximation in the case of Schrodinger operators |13| and coincides with the small-r 
asymptotics of the GUE-form factor ( |2.21| ). If, however, the quantum Hamiltonian 
is invariant under time reversal so that we have to choose g = 2, the diagonal 
approximation becomes K^^ir^Tu) ~ 2r. According to ( |2.17D Kramers' degeneracy 
then forces us to compare the random matrix form factor with the modified semiclassical 
result 

Kf aa { T] T H ) = l -Kf a9 (I; T H ) ~ ~r , (5.1) 

which now coincides with the small-r asymptotics of the GSE-form factor ( 2.19|) . In 
both cases this is exactly the behaviour that is predicted by the conjecture of Bohigas, 
Giannoni and Schmit as stated in ( [2.20|) and (|2.18|) . 
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Appendix A. Classical and quantum mechanical time reversal 

In this appendix we are going to discuss the relation between classical and quantum 
mechanical time reversal for systems with spin 1/2. We restrict our discussion to the 
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conventional time reversal operation that leaves the position coordinates unchanged 
and reverses momentum and spin coordinates as well as the time t. The discussion 
of generalized time reversal operators that combine conventional time reversal with 
another unitary symmetry operation is analogous, see 0] for examples. To be specific 
we consider as a quantum Hamiltonian for a particle of mass m, charge e and spin 1/2 
either a Dirac-Hamiltonian 

H D = col f y V - ^ A i x )j + fi m( ? + eip(x) , (A.l) 

or a (generalized) Pauli-Hamiltonian 

H P = H s t 2 + h(rC \^\7,x J . (A.2) 



x i 

In the relativistic case (|A.1|) the Dirac algebra is realized by the 4x4 matrices 




« = ( : ; i . = i i _° a 1 ' (A,3) 

where cr denotes the vector of Pauli matrices and I2 is a 2 x 2 unit matrix. The non- 
relativistic Hamiltonian ( |A.2j ) is composed of a Schrodinger operator H$ and a coupling 
term of spin to the translational degrees of freedom. The latter has to be understood 
as the quantization of some M. 3 - valued function C(p,x) on phase space. For example, 
this can be a magnetic field, i.e., Cb(p, x) = B(x), or a spin-orbit coupling term 
C so (p, x) = 4m2 ].2 lxl dV J^ l) (x x p). 

For systems with spin 1/2 the operator of time reversal is given by 

f := e'i^ K = i o y k , (A.4) 

where K is the operator of complex conjugation in position representation, see In 
the relativistic case, where T has to act on four component spinors, ( A.4 ) shall mean a 
block diagonal 4x4 matrix with two copies of ( |A.4|) in the diagonal blocks. A quantum 
system with Hamiltonian H to be time reversal invariant requires THT^ 1 = H. Thus, 
in the case of a Dirac-Hamiltonian ( |A. 1|) 

TH D f- 1 = c{-a) ( --V - -A(x)) + (3mc 2 + ep(x) (A.5) 



1 c 

shows that conventional time reversal invariance is equivalent to the absence of magnetic 
forces. For the Pauli-Hamiltonian ( |A.2| ) to commute with T we first need the Schrodinger 
operator Hg to be time reversal invariant, i.e., KHsK = Hs- In addition, the condition 

faC (j V, x\ f- 1 = -<tC x \ = aC (j V, x^j (A.6) 

has to be met, i.e., the coupling term C(p,x) must be an odd function of momentum 
p. This requirement is fulfilled by the spin-orbit coupling term C so , but is violated by 
the coupling Cb to an external magnetic field. In both the relativistic and the non- 
relativistic situation, however, even the presence of a magnetic field might allow for the 
existence of an anti-unitary operator representing a generalized time reversal symmetry 
that commutes with the Hamiltonian, see IBL [2(1 . 
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We now want to investigate the implications of a quantum mechanical time reversal 
invariance for the semiclassical analysis of the form factor. The first obvious consequence 
is a time reversal invariance of the classical translational dynamics. But furthermore, 
also the spin dynamics, governed by the spin transport equation ( |4.6| ), exhibits a certain 
kind of symmetry under time reversal. In order to discuss the latter, one has to study 



the behaviour of the matrix M entering (|4.6|) under x i— ► x, p i— > —p. For a Pauli- 



Hamiltonian M is given by crC(p,x), where C(p,x) is defined as above, and for a 
Dirac-Hamiltonian with no magnetic field one finds M = g(\p\)x X p, see |15[ for 
details. In each case quantum mechanical time reversal invariance therefore implies 
M(—p, x) = —M(p, x). Substituting now tt—>—t and p i— > — p in ( [4.6|) one obtains 



- d\-p, x, -t) + i M(^ H \-p, x)) d(-p, x, -t) = . (A.7) 

Since the translational dynamics are time reversal invariant, M being odd in the 
momentum variable leads to 

d(-p, x, -t) + i M{¥ H (p, x)) d(-p, x, -t) = , (A.8) 

and therefore d(—p, x, — T 7 ) = d(p, x, T 7 ). Moreover, the fact that d(p, x, T 7 ) G SU(2) is 
a holonomy implies d(—p, x, — T 7 ) = [d(—p, x, T 7 )] _1 so that finally tr d(—p, x, — T 7 ) = 
trd(p, £E, T 7 ). Altogether the above considerations confirm that the presence of a 
quantum mechanical time reversal invariance implies that a primitive periodic orbit 
7 and its time reversed partner share identical amplitudes A,,^. 

Appendix B. Equidistribution of long periodic orbits 

In this appendix we want to show how the periodic-orbit representation ( [4. 4ft of Liouville 
measure can be obtained from equidistribution properties of periodic orbits. The basic 



reference for the following is [25]. Our assumptions on the flow $^ are as stated in 
section f|. 

In order to proceed further we first have to introduce some notation. Let B denote 
the set of all "^-invariant Borel probability measures on Qg. Any \i 6 B can be 
associated a metric entropy h^. Then for any Holder-continuous observable / G C a (flE), 
with some a > 0, the topological pressure is defined as 



P(f) := sup |^ + f dfi ; fi G Bj . (B.l) 

The supremum is attained for a unique measure fif, which is called equilibrium measure 
for the observable /. The periodic orbits of the flow $^ are then equidistributed with 
respect to /// in the following sense, 

a(p, x) d/i/(p, x) = km ^ — f , (B.2) 

for every a G C(Qe)- The averages over periodic orbits are defined as in (|4.5|). 
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In a next step one has to identify Liouville measure as the equilibrium measure of 
some observable /. To this end we consider the tangential map D& H restricted to the 
unstable subbundle E u of the tangent bundle TQe and define 

f(p,x) := -Al og det D¥ H (p,x)\ Eut=Q . (B.3) 
A direct calculation then yields 

^ d-l 

p = -—J2u 1 , j so that e T - /7 =p 7 , (B.4) 

7 3=1 

compare (|4.3| ). Furthermore, the equilibrium measure associated with the observable 
( |B.3|) is called Sinai- Ruelle- Bow en measure hsrb, for which it is known that for every 

a G c(n E ) 

lim i / a(& H (p,x)) dt= [ a(p',x') dfx S RB(p',x') (B.5) 

T^oo 1 J Q J Qe 

holds for a set of initial conditions (p, x) G with positive Liouville measure, see 
| 29fl for more information. Since is supposed to be ergodic with respect to one 
concludes that Hsrb — Me- 

What is still lacking is an asymptotic estimate of the denominator on the right-hand 
side of ( |B.2|) . In order to obtain this we first have to introduce a regularization in that 
we multiply the observable (|B.3| ) by some factor (3 <1. Then we appeal to the relation 

? ^"iwi''" 1 "^ 11 ■ T "°°' (R6) 

T-e<T 7 <T+e 



see ||25|| . Since P(f) = 0, in the limit (3 — > 1 one obtains for k G N 

T*p 1 ~2eT k -\ T^oo. (B.7) 

7 

T-e<T-,<T+£ 

Replacing now s by T, followed by the rescaling 2T i— > T, this implies 

T * Pi ~ 2^ fc T k , T -> oo . (B.8) 

7 

T 7 <T 

Introducing then the relations ( p.4j) and (p.8|) in (|B.2|) finally yields the periodic-orbit 
representation (f4.4|) of Liouville measure. 
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